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Abstract 

Principles of self-organization play an increasingly central role in models of 
human activity. Notably, individual human displacements exhibit strongly recur- 
rent patterns that are characterized by scaling laws and can be mechanistically 
modelled as self-attracting walks. Recurrence is not, however, unique to human 
displacements. Here we report that the mobility patterns of wild capuchin mon- 
keys are not random walks and exhibit recurrence properties similar to those of 
cell phone users, suggesting spatial cognition mechanisms shared with humans. 
We also show that the highly uneven visitation patterns within monkey home 
ranges are not entirely self-generated but are forced by spatio-temporal habitat 
heterogeneities. If models of human mobility are to become useful tools for pre- 
dictive purposes, they will need to consider the interaction between memory and 
environmental heterogeneities. 

Keywords: movement ecology, home range, memory, human mobility, 
capuchin monkeys, scaling laws. 
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1 Introduction 



Individual human displacements have a strong impact on many collective social phe- 
nomena, such as the spread of epidemics [1, 2, 3] or cultural traits [4]. The availability 
of quantitative mobility data has increased in recent years through the widespread use 
of global positioning systems [5] and the ability to track cell phones [6, 7, 8]. These 
new data show that human displacements do not follow memoryless processes, like 
the well-known Markovian random walk [9], rather, they exhibit ultra-slow diffusion 
and unusual, long lasting recurrence properties due to the tendency of individuals to 
frequently revisit a small number of familiar locations [7, 8] . 

Whether unifying principles that govern the movements of humans [8] and, more 
generally, of living organisms [10] exist is a hotly debated issue. As with the distribution 
of city sizes [11] or the dynamics of individual tasks [12], scaling laws prevail in human 
mobihty data [5, 8, 13]. In physical systems, scaling laws are often the outcome of self- 
organization principles. Self- attracting or reinforced walks [14, 15] are non-Mar kovian 
stochastic walks that tend to revisit with higher probability locations visited in the 
past. Recently, a new reinforced walk model with preferential return to previously 
visited sites, in analogy with preferential attachment rules used in network science [16], 
could reproduce many empirical human visitation patterns [8]. However, it is also 
important to quantify the dependence of individual displacements on more complex 
external factors, for instance, transportation [1] and social networks [17], or the spatial 
distribution of facilities or resources [18]. 

In this contribution, we show that the scaling laws characterizing recurrence in hu- 
man movements are similar to those exhibited by some foraging animals, which suggests 
that these patterns are not unique to humans and may be generated by a more general- 
ized set of cognitive mechanisms. We present evidence that, like humans, other animals 
also live in heterogeneous habitats where some areas are more valuable to them than 
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others. Recurrent patterns of movement, in particular home range behavior, actually 
characterize the ranging of a very large number of animal species [19, 20]. Here, we use 
radio telemetry data to show that the mobility and visitation patterns of capuchin mon- 
keys ( Cebus capucinus) foraging for fruit in a tropical forest are not only comparable to 
human cell phone users qualitatively, but exhibit scaling laws that are strikingly similar 
in both the short and long time regimes. In addition, monkeys movements are not en- 
tirely self-organized but also strongly driven by spatio-temporal variations in resource 
distribution. 

2 Results 

The data for our analyses come from a study on Barro Colorado Island, Panama, in 
which the movements of four radio-collared capuchin monkeys belonging to different 
social groups were tracked using an automated telemetry system over the course of a 
six month period, from November 2004 to April 2005 (see [21] for details). The locations 
for each individual were estimated every A^o = 10 minutes and discretized into square 
cells of size 50x50 m, a scale that corresponds roughly to the measurement error and was 
taken as the spatial resolution. Unless indicated, we considered all available positions 
(night and day) in the analyses. 

2.1 Basic movement properties 

Short time-scale properties of monkey movement paths followed similar scaling forms 
as those of human cell phone users Ref. [8] (see Figure la,b). For example, the dis- 
tributions of the displacement Ar between consecutive telemetry fixes are good fits to 
a power-law, (Ar)~^~" with a ~ 1.7 in the range 50 m < Ar <350 m, for all four 
individuals (sample sizes=1619, 1998, 3072 and 2097). This exponent, as well as (3 and 
rj below, were obtained by using a Bayesian Markov chain Monte Carlo simulation with 
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Figure 1: Comparison of monkey and human movement properties. Data for each of 
four radio-collared capuchin monkeys plotted separately. Monkey power-laws exponents 
(subscript "m") of the frequency distributions are estimated from pooled data using 
Bayesian Markov Chain Monte Carlo Estimation. These values are practically indistin- 
guishable from the MLEs, which are not shown. Human power-law exponents from [8] 
(subscript "h") are given in each plot for comparison, a) Displacements between consec- 
utive telemetry fixes (Ar in meters) showed a much steeper decay, P(Ar) ~ (Ar)^^^° 
with larger a, than in humans but b) waiting times (1 unit=10 min) showed a power- 
law exponent close to humans, c) RMS displacement ((r(t) — r(0))^)^/^ (in meters) 
rose quickly during the first day but d) very slowly thereafter, with values always much 
smaller than the size of Barro Colorado Island (diameter~ 5km). e) The number of 
distinct 50x50 m habitat cells visited, S{t), as a function of time, f) Like for humans, 
the probability that a cell received n visits within the 6-month time interval, P{n), 
decays as a power-law, but with a smaller exponent. Estimation cutoffs for Ar, At and 
n were set respectively at 350 m, 12 h, and 35 visits. In all cases the bounded power-law 

model was a very good fit to the data (i?^ > 0.98) and the frequency distribution with 
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MLE had an Akaike weight > 0.9999 when compared to the most likely exponential 
distribution. The sizes of the 95% confidence intervals of the exponent values were 
always lower than 0.20. 



uniform priors and the Metropolis-Hastings algorithm. For each measure we evaluated 
the likelihood of a given value of the exponent of the bounded power-law given the data. 
The exponent value is the average obtained from the posterior support distribution. We 
also calculated the maximum likelihood estimate (MLE), obtaining practically the same 
value. 

Similarly, the waiting time At, defined as the number of consecutive fixes found in 
the same 50x50 m cell multiplied by the temporal resolution Ato, follows a power- law 
distribution, P{At) ~ {At)~^~^ with /3 = 0.8, the same value as reported for humans 
[8, 13] and other primates [22]. The exponents a < 2 and /3 < 1 play central roles in 
continuous time random walk (CTRW) models of anomalous transport [13, 23] and of 
biological Levy flights [22, 24, 25]. 

Like humans, monkeys trajectories are however incompatible with a random Levy 
walk description at large time-scales, as they exhibit strongly recurrent patterns of 
movement. For example, the root mean squared displacement (RMSD) virtually asymp- 
totes within one day, rather than following a classical power-law growth with time (Fig. 
lc,d). This arrested or ultra-slow diffusion is a manifestation of home range behavior 
[19, 20]. Likewise, the average number S{t) of distinct cells visited during a time in- 
terval t (Fig. le) increases through time as with fj, smaller than the waiting time 
exponent /3, which is what would have been expected for a CTRW [8, 23]. One finds 
fj, — 0.6, remarkably close to what measured by Song et al. for humans [8]. Another 
intriguing aspect of recurrence is the highly uneven visitation pattern among cells, see 
Fig. If. Within the six-month data collection period, the probability P{n) of finding a 
cell that has received n visits by the same individual {i.e., has been entered n times) 
cannot be described by a Poisson distribution or a bell-curve centered around a char- 
acteristic value. Instead, it can be fitted by a power-law, P{n) ~ n~'^ with 77 = 1.0 
{rj — 1.8 in humans [8]). Whereas most sites were visited once or twice, it is not rare 
to find popular activity "hotspots" or "hubs", with 30 or more visits. In contrast, a 
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random walker bounded in a closed domain would visit all sites more or less equally 
frequently. For all the properties described above, the four individual monkeys showed 
very similar exponent estimates. The number of distinct cells visited by the animals 
during the six-month period was 568, 488, 606 and 531, respectively. 

2.2 Recurrence and environmental forcings 

Like humans, many animals have sophisticated memory skills [26] and home ranges can 
emerge from frequent returns to previously visited locations [27]. In principle, home 
ranges could even be completely self-organized, arising as a result of the walker's history. 
Random walks biased towards the center of mass of all previously visited sites, or such 
that, at each step, there is a finite probability of returning to a randomly selected 
site visited in the past, exhibit arrested or very slow diffusion [20, 28]. In the context 
of human mobility, the results of Figure 1 have been reproduced assuming that the 
probability of choosing a given site is proportional to the number of previous visits to 
that site (''preferential attachment"), and assuming that the probability of taking a 
random step (to a unvisited site) decays algebraically with the number of sites already 
visited, S{t) [8]. An advantage of this latter approach is that it builds up uneven visits 
and power-laws for the distribution P{n). 

What these models and those belonging to a much broader class based on preferential 
attachment mechanisms ignore is environmental heterogeneity. A mobile agent may be 
inclined to visit a place not because it is familiar but because of some intrinsic quality of 
the location, {e.g., food content). Heterogeneity is implied in the relationship between 
the quality of habitat cells and the probability of cell revisit. Capuchin monkeys home 
ranges typically contain millions of tree stems which belong to more than 200 species 
[29], vary in size by 2 orders of magnitude [30] and only episodically produce fruit [31]. 
Consequently, the cells in our analyses varied widely in the content of monkeys' primary 
food, ripe fruit. During the day (from 04:00 to 20:00 hrs), animals are awake and active 
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[21], and their primary activities are feeding and foraging with relatively rare periods 
of prolonged resting [29]. During the day, a waiting time of At = Ato = 10 min (one 
fix, the minimum observable value) is assumed to correspond to a "transit" through a 
cell. Given a series of consecutive fixes in a same cell, the time elapsed between the first 
and the last fixes. At — Ato, can be taken as a proxy of the ripe fruit content within 
the corresponding cell. An individual will tend to spend more time feeding in a tree 
that has more fruit and we consider At = Ato as the threshold above which one can 
tell that a feeding event occurred. Figure 2 shows a positive correlation between the 
number n of day-time visits received by a cell and the average duration of a visit at 
that cell, {At — Ato)n (the average being taken over all visits among cells with n visits). 
The latter quantity can be approximated by a scaling law: 

{At - Ato)n ~ n\ (1) 

with 7 f« 0.3. Although other functional relations may fit the data, the crucial point 
is that large n cells have, on average, longer waiting times and, most likely, higher 
food content. A similar, and even clearer, relationship is observed if all night and day 
positions are taken into account (7 ^ 0.4 in this case). Sleeping trees also represent 
resources, as not all trees are suitable for sleep at night. 

Similarly, the time elapsed between two consecutive day-time visits to a site depends 
on the results of the initial visit. We denote MFRTa* as the mean first return time 
to a cell where the waiting time was At at the last visit. As shown by Figure 2b, the 
MFRT tends to decrease with increasing patch quality (At). For At < 100 min, it 
approximately follows the scaling law: 

MFRT At - (At)-^ (2) 

with 6 ~ 0.3. Note that the unaveraged first return time is a variable known to fluctuate 
widely in simple models like the random walk [32]. Similarly here, error bars in Figure 
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Figure 2: a) Average visit duration at a cell that has otherwise received n visits during 
a time window corresponding to 5" = 200 distinct visited sites (about 2 weeks period). 
Each point represents an average over sites with n visits in one time window, over 12 
different time windows and the 4 different monlceys. The curve obtained from all fixes 
("all data") has been upscaled by a factor 2 for clarity, b) Mean first return time to a 
cell as a function of the time spent on that cell during the last visit. Averages obtained 
from less than 5 observations were discarded. The solid lines are least-square fits and 
have slopes: 0.41 and 0.27 in a) (all data and day time data, respectively); —0.39 and 
—0.29 in b) (all data and day time data, resp.). In all cases, the probability that the 
variables are uncorrelated {S — ^ — 0) is p < 10~^. 
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Figure 3: a) Day-time (4:00-20:00hrs) positions occupied by the same individual monkey 
during successive, non-overlapping periods each containing S = 200 visited sites. The 
black dots represent the 10 most visited cells during the day, for each time period. 
Chronological order goes from left to right and top to bottom. Hotspot probability Uhs 
(b) and mean number of visits to a hotspot N^s (c), as a function of the separation m 
between time windows (see text). 
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2b increase with At, which may be due to the fewer available observations at larger At 
(Fig. lb). Averages obtained from less than 5 observations were discarded. 

These results show evidence that, in addition to being non-random, movement 
choices are driven by resource availability, with larger food patches visited more of- 
ten. Eq. (2) also holds if night and day positions are aggregated (7 0.4). Hence, 
a new waiting time is not randomly drawn at each cell visit, as assumed in CTRW 
models. Rather, the situation here is analogous to systems with quenched site disorder 
[32]. In trap models, though, the frequency of visits to a site by a memoryless random 
walker is independent of the waiting time {'-f — 5 — 0). 

Unlike supermarkets, gas stations or other human facilities, resources in natural 
ecosystems tend to be ephemeral. Trees produce fruit transiently and often asyn- 
chronously, with the fruiting periods of individual trees lasting from a few weeks to 
a couple of months [31]. Capuchin monkeys hotspots, in particular, are not permanent. 
Figure 3a shows the positions of the 5% most visited sites by a single individual during 
the day in successive non-overlapping time periods (of about 2 weeks) with S — 200 
visited sites each. These cells are rarely the same from one time window to the other. 
The fluctuating cloud formed by the visited cells suggests a flexible use of space, despite 
the stability of the overall home range. Defining hotspots as the most frequently visited 
10% of sites, the probability Uhsim) that a cell which is a hotspot during a time window 
i is also a hotspot during time window i + m decays rapidly with m (Figure 3b). The 
average number of visits Nhs{m) received by a cell in time window i + m, given that 
this cell was a hotspot in time window i, also decays with m (Figure 3c). Very similar 
curves are obtained if night and day data are aggregated. 
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3 Conclusions 



We have shown that the individual movement patterns of humans are similar to those 
of capuchin monkeys. Most notably these patterns exhibit ultra-slow diffusion (the 
presence of a home range) and history dependent recurrence properties obeying scaling 
laws. In monkeys, we have presented evidence that these nonrandom walks are caused 
by the use of memory, a behavior which is not unique to humans nor evolutionary novel. 
These results suggest, but do not prove, that similar cognitive mechanisms may govern 
the movements of many animals. Testing such hypothesis would require data analysis 
for a variety of species. 

Our analyses also indicate that monkey movements are driven by environmental 
heterogeneities. Monkey ranging patterns are thus not "self-quenched" into a routine 
emerging from initially random movements that would be re-inforced and dominated 
at large time by frequent revisits to a small number of known locations, a mechanism 
proposed for human mobility [8]. Similarly, it is likely that human movements are not 
entirely self-organized, as this would imply that the locations of activity hotspots would 
be uncorrelated with environmental factors. In the real world, many hotspots would 
clearly be shared because of some intrinsic property of the location: e.g., schools, transit 
hubs or office buildings. A difference is that many human resources last over long time 
scales, like homes and workplaces [6]. Yet, others are transient, like restaurants, trendy 
night-clubs or fashionable shops, or seasonal, like swimming pools and ski-resorts. A 
place can be also abandoned for another not because it is depleted but because the 
appetite for that resource has been temporarily sated. 

Despite the similarities between human and monkey movements, there are important 
differences. In humans, the home range size fluctuates widely from one individual to 
another, ranging between 1 km and 1000 km approximately [6, 8]. In contrast, the 
nearly identical home range sizes found for the four monkeys (Fig.lc,d) suggest a narrow 
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size distribution, probably due to comparable habitats, energetic needs and locomotor 
capacities. This is consistent with the fact that the home range of a capuchin group 
overlaps with a relatively small number of other home ranges [21], which probably 
limits intergroup interactions to nearest neighbors or next-nearest neighbors. This 
property is not shared by all animal species, though. Territorial prides of Serengeti 
lions occasionally make long-range contacts with other prides as well as with nomadic 
individuals, producing a denser network of contacts with the "small world" property 
that characterizes many human social networks [33]. 

Similar to random diffusion models of animal movements with intermittent cen- 
tral attraction [34] , several models of human mobility with preferential return to home 
have been introduced recently, showing that recurrence can have dramatic effects on 
spreading processes in large populations [1, 2, 3]. These latter models typically as- 
sume a known set of locations that can be visited and Markovian, individual-dependent 
transition rates between locations. Developing parallel modeling frameworks that ad- 
equately represent both self-organization and environmental forcing will be critical to 
the success of human movement models for purposes such as of controlling the spread 
of infectious diseases. These alternate models should incorporate cognitive mechanisms 
that many vertebrates use, e.g. spatial representation mechanisms (cognitive maps or 
travel cost discounting) and temporal mechanisms (episodic memory) [35]. How scal- 
ing laws emerge from the interplay between memory and landscape features remains 
elusive. 
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